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We study the probability distribution function (PDF) of the smallest eigenvalue of 
Laguerre-Wishart matrices W = where A is a random M x N (M > N) ma¬ 

trix, with complex Gaussian independent entries. We compute this PDF in terms of 
semi-classical orthogonal polynomials, which are deformations of Laguerre polynomials. 
By analyzing these polynomials, and their associated recurrence relations, in the limit of 
large N, large M with M/N ^ 1 - i.e. for quasi-square large matrices X - we show that 
this PDF, in the hard edge limit, can be expressed in terms of the solution of a Painleve 
III equation, as found by Tracy and Widom, using Fredholm operators techniques. Fur¬ 
thermore, our method allows us to compute explicitly the first 1 /N corrections to this 
limiting distribution at the hard edge. Our computations confirm a recent conjecture by 
Edelman, Guionnet and Peche. We also study the soft edge limit, when M —A ~ 0 {N), 
for which we conjecture the form of the first correction to the limiting distribution of 
the smallest eigenvalue. 


1. Introduction 


The study of extreme eigenvalue statistics in Random Matrix Theory (RMT) has 
attracted much attention during the last twenty years. In particular, the Tracy- 
Widom (TW) distributions |1I2| describing the largest eigenvalue Amax (as well 
as the smallest one, Amin) in the classical Gaussian ensembles, orthogonal (GOE, 
/3 = 1), unitary (GUE, /? = 2) and symplectic (GSE, /3 = 4) - where /3 is the 
Dyson index - have become cornerstones of the theory of extreme value statistics 
of strongly correlated variables. Quite remarkably, it was shown that the TW distri¬ 
butions, denoted by rp in the following, appear in a wide variety of problems m , a 
priori not directly related to RMT, ranging from the longest increasing sequence of 
random permutations of integers jj], stochastic growth and related directed poly¬ 
mer models in the Kardar-Parisi-Zhang |5] universality class |6l7l8l9linill] and 
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sequence alignment problems m to non-intersecting interfaces and Brownian mo¬ 
tions | 13|14|15| as well as in finance m- While the TW distribution describes the 
typical fluctuations of Amax and Amin, a large body of work has also been devoted 
to the study of large deviations of extreme eigenvalues in Gaussian ensembles HU. 


soft edge hard edge 




N N 

Figure 1. Plot of the Marchenko-Pastur distribution for c < 1 (left) and for c = 1 (right). 


Another interesting ensemble of random matrices, which we focus on in this 
paper, is the so-called Wishart-Laguerre ensemble - here we focus on the case of 
complex matrices (/3 = 2). Let X he a, M x N rectangular matrix with i.i.d. complex 
Gaussian entries and M — N = a > 0. The Wishart-Laguerre matrix W is defined 
as W = X which is thus a N x N Hermitian matrix, having N real and positive 
eigenvalues Ai, A 2 , • • • , \n- The joint probability density function (PDF) of these 
N eigenvalues is given by HHnn] 


N 


N 


N 


Pjoi„t(Ai, A 2 , ■■■ = — n(A. - A,)" exp ( - ^ A, ) n A“, (1.1) 


i<3 


where Ai, A 2 ... Xn are positive and where Zjq is a normalization constant, depend¬ 
ing on a. These Wishart-Laguerre matrices play an important role in statistics, in 
particular in principal component analysis, where the matrix W is a covariance ma¬ 
trix. Hence in this case both M and N, and thus a = M — N, are positive integers. 
However, this joint PDF in Eq. (1.1) is well defined for any real value of a > 0, 
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which, for some non-integer value of a, may have some physical applications. An 
interesting example is the case of non-intersecting Brownian excursions |2()| . i.e., N 
non-colliding positive Brownian paths Xi{t) > 0 on the unit time interval t S [0,1] 
constrained to start and end at the origin Xi{0) = Xi{l) = 0. The joint PDF of the 
positions of the N walkers, at a given time t, can indeed be written as |2T] 

^ N N ( ^ 

Pioint{xi, ■■■ ,XN]t) = y-- xf W(.xl - xjf exp I - ^ 

^ i—1 i<j \ i—1 



where cr^(t) = 2t{l — t) and ZN{t) is a normalization constant. From this expres¬ 
sion (1.21 we obtain that the scaled variables xf/a'^{t) behave statistically like the 
eigenvalues of random matrices from the Wishart-Laguerre ensemble with a 
non-integer parameter a = 1/2. Hence it is physically relevant to study the joint 
distribution in Eq. 0 for any real a > 0. 

A first important characteristic associated to this ensemble (1.11 is the mean 
density of eigenvalues, p{X,N), defined by 


1 


N 


^ 2 = 1 


(1.3) 


where the overline denotes an average over the different realizations of the random 
variables A^’s according to the joint PDF in Eq. 0. In the large N limit, it is 
well known that p{X,N) is given by the Marchenko-Pastur (MP) distribution (see 
Fig. 


p{X,N) 


1 

—^ Tr 

N—¥oo N 


Pmp 


Pmp (x) = 




27ra; 


(1.4) 


(1.5) 


where x± = ± 1)^ are the right and left edges of the support, with c = 

N/M < 1. The case where c < 1 corresponds to the case where a ^ 0{N). Here, 
we mainly focus on the case c = 1 which corresponds instead to the case where a 
is finite, while both N and M are large. In this case, which we will mainly focus 
on in this paper, the MP distribution takes the particular form (see the left panel 
of Fig. 

, , 1 U-x 

Pmp \x) = 


27r 


( 1 . 6 ) 


At the right edge, near = 4, pmp (x) vanishes as a square-root, pmp (x) oc — x 
and therefore the fluctuations near this soft edge are governed, for large N, by 
the Airy kernel. In particular, the distribution of the largest eigenvalue Amax) 
appropriately shifted and scaled, converges, when N —>■ oo, to the TW distribution 
p 2 mentioned above (the same as for GUE) |6I22| . It is now well known that this 
distribution can be expressed in terms of a special solution of a Painleve H equation. 
While this connection to Painleve transcendents was initially obtained by Tracy and 
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Widom using Fredholm operator techniques, Nadal and Majumdar |2H| provided, 
more recently, a derivation of this resnlt (for (5 = 2) nsing semi-classical orthogonal 
polynomials (OPs), see also Ref. [21]. This method is at the heart of the present 
paper. 

On the other hand, at the left edge, X- = 0, pmp {x) has a sqnare-root singu¬ 
larity, pmp {x) oc \fx (see the left panel of Fig. [^. What about the fluctuations 
of the smallest eigenvalue Amin hr this case? One can estimate the typical scale 
of Aminj for large N, by considering that there is typically one eigenvalue in the 
interval [0, Amin], he.. 


N 


p{X,N) dA~(P(l) , 


(1.7) 


which implies that Amin ~ 0{1/N). Within this scale 0{1/N), the fluctnations are 
governed by the Bessel kernel |25I26| . Fnrthermore, it has been shown that the 
distribution of A^Amin converges to a limiting form which (i) is different from the 


standard TW distribution J-<i and depends continnonsly on the exponent a in (1.11 


and (ii) can be written in terms of a special solution of a Painleve III eqnation m- 
If one introduces FN{t) = Pr(Amin > t) then, one has indeed 

f{u) 

^ N \ — \ = = tlXP [ / 

N—>-oo 


lim iW( ^ ) = F^{x), F^{x) = exp ( / 

\N / \Jo 

where f{x) is the unique solution of a Painleve III |27| : 


-dit 


( 1 . 8 ) 


{xfy+^r{i+r){xf-f) = {afr 


satisfying 


fix) ~ - 


j,a+l 


r(a-k l)r(a-f 2) 


as a: —?► 0 . 


(1.9) 


( 1 . 10 ) 


This result (1.8 1 was shown by Tracy and Widom using Fredholm operator tech¬ 


niques m- Note that for integer values of a, the limiting distribution F^dx) can 
be written as an a x a determinant whose entries are expressed in terms of Bessel 
fnnctions (see Eq. (C.l I below) - a resnlt which can be obtained by clever manip¬ 
ulations of determinants |2H] • In particular for a = 0 the result is extremely simple 
as Fpf{x) = exp {—Nx) for all N, implying /(x) = —x, which is obvionsly solution 


of Eq. (1.9) with the boundary condition (1.10). 


The limiting distribution of An 


N —>■ oo, is thns well known (1.8) 


for complex Wishart matrices, in the limit 
we refer the reader to Ref. [22] to a recent 
work on the smallest eigenvalne for real Wishart matrices in the hard edge limit. 
What abont the finite N corrections to this asymptotic form? Such a question 
is quite natural for practical applications of extreme value statistics (EVS), where 
one always deals with finite samples - here matrices of finite size. This issne was 
recently revisited for EVS of independent and identically distributed random vari¬ 
ables nsing a renormalization gronp approach m- For EVS of strongly correlated 
variables, there are actually few cases where these corrections have been worked ont 
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explicitly, including random walks m, the largest eigenvalue of random matrices 
belonging to various ensembles non-intersecting Brownian motions 

pH] or (Poissonized) random matchings [ 37 ]. For real Wishart matrices, the first 
corrections to the limiting distribution of the smallest eigenvalue in the soft edge 
limit were studied in Ref. |85| where it was shown that corrections to the limiting 
distribution of Amin and Amax are quite different, although the limiting distributions 
for both observables are actually the same, namely the TW distribution for GOE, 


Fi- What about the corrections to the limiting distribution in Eq. (1.8 1 of Amin 
for complex Wishart matrices in the hard edge limit? This question was recently 
raised by Edelman, Guionnet et Peche |3S] in their study of finite size covariance 
matrices with non-Gaussian entries. Based on the large N expansion of the exact 
formulas obtained in Ref. jSS] for small integer values of a, they conjectured the 
following form of the first 1 /TV-correction 

=E^(x) + ^xF,^{x) + o(^) . (1.11) 


Note that this first 1/iV correction in Eq. (1.11) can be interpreted as a correction 
to the width, i.e.. 




( 1 . 12 ) 


It is interesting to notice m that for most of the cases which have been studied 
in RMT |82l88l84| . it was actually found that the first order correction to the 
limiting distribution of extreme eigenvalue actually corresponds to a correction of 


the scaling variable, as in Eq. (1.121. One exception concerns the smallest eigenvalue 


of real Wishart matrices in the soft-edge limit, where the first correction has a more 
complicated structure |85| . 

The main goal of this paper is to provide an explicit computation of this first 
correction in the hard edge limit and we will show that it has indeed the conjec¬ 


tured form given above in Eq. (1.11). To perform this computation, we will use a 
method relying on semi-classical OPs, in the spirit of Refs. [ 33 ] and [ 33133 ] . As we 
will see, our method does not only allow us to compute explicitly the first 1/N cor¬ 
rections but provides also a rather straightforward derivation of the expression for 
the limiting distribution F’oo(x) in terms of the solution of a Painleve III equation, 
without using Fredholm operators theory but relying instead only on the recurrence 
relations associated to the (semi-classical) OPs system. Finally, we will also study 
the first finite N corrections to the limiting distribution of Amin at the soft edge. 

Note that after the results obtained in the present paper were presented in a 


conference m, another independent proof of the conjecture in Eq. ( jl.llj ) was 
achieved in Ref. lu, using operator theoretic techniques. More recently, yet an¬ 
other independent proof of this conjecture was given in Ref. |42| . 



















6 


2. Summary of main results and outline of the paper 


The distribution of the smallest eigenvalue Amin = min A^ is given by 

l<i<N 

/ oo poo poo 

dXi dX2- dXNPjoint{Xi,X2,- ■ ■ ,Xn) ■ 

( 2 . 1 ) 

In this paper, we will compute using semi-classical OPs {7rfc(A)}fcgN which 

are polynomials of the variable A while t and a are parameters (for the sake of 
clarity in the notations, this dependence is omitted here): 


(TTfelTTfe/) = /“e ^A“7rfc(A)7rfe/(A)dA = 


( 2 . 2 ) 


T^kiX) — A* -|- CfcA* ^ -f ... 

The cumulative distribution Fj^(t) can be expressed, using standard manipulations, 
in terms of the norms hk’s as 

N-l 


N\ T-r 

FN{t )=^ n ■ 


Zn 


(2.3) 


fc=0 


As we will see, the norms h^’s can be computed from the three term recurrence 
relation satisfied by the OPs 

Att^ = -t- SkTTk F Rk'^k—i , (^-4) 

from which we deduce the following important relations: 

Rk = (2.5) 

f^k-1 

Sk = —tdt log /^fc + 2A: + a + 1 , (2.6) 

k-l 

Ck = -Y.S,. (2.7) 

i=0 

Note that the starting point of our analysis is very similar to the one of Basor 
and Chen in Ref. [31] but the analysis of the recursion relations is different. In 
particular, we do not make use of ladder operators techniques, which are heavily 
used in Ref. [21]. In addition, we provide an asymptotic analysis of this OP system 
(2.21 for large N, beyond the leading order. 

In section]^ we will study the variables Sk, Rk, hk and (^k- In particular, we 
will show that Sk and Rk satisfy a coupled set of equations, named Schlesinger 
equations in the literature on OPs 

Sk — Rk+i + Rk = tdtSk, 

( 2 . 8 ) 


2 — Sk+i + Sk — t 


dtRk+i 

Rk+l 


This system of equations, together with the initial condition given in Eq. (3.331 
below, determines uniquely the values of Rk and Sk for all values of fc, as they 
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can be computed by induction. It is however quite difficult to analyze the large N 


behavior of Rn and Sn using only this set of equations (2.8). To circumvent this 


difficulty, it is customary to use another set of relations, called the Laguerre-Freud 
equations, which we derive here using the method based on Turan determinants, 
as developed in Ref. (43] . They read 

Rk +2 ~ Rk = 5'fe+i(2fc + 4 + a + t — Sk+i) — Sk{‘ik + a + t — Sk) — 2t, 

Sk+i{Sk+i ~ t) = Rfe+i(2fc + 1 + a +1 — Sk+i — Sk) ^ ^ 

~ Rfe+2(2fc + 5 + a +1 — Sk+2 ~ Sk+i)- 

Some details of this derivation have been relegated in [Appendix A[ 

In section]^ by manipulating these two sets of equations (2.8 1 and (2.9), we show 
that Fiy{t), for finite N, is related to a special solution of a Painleve V equation [see 
Eqs. ( |4.3| 4.15)], thus recovering a previous result of Tracy and Widom. Section]^ 
is devoted to the hard edge scaling limit: 

N^oo,t^0,x = NtG fixed . (2-10) 

It is known that the limiting distribution Fao{x) can be obtained by analyzing the 


large N limit of this Painleve V equation in Eqs. (4.3 4.15), obtained for finite N, 


leading to (1.8) and (1.9). As we show in section 5, the finite N corrections are 


then easily obtained from the Schlesinger equation (2.8), from which we obtain Eq. 
(1.11). Finally, section|^is devoted to the soft edge scaling limit, when a ~ 0{N), 
for N large. 


3. Semi-classical Orthogonal Polynomials 

To study this OP system (2.2) it is useful to introduce a deformation parameter 
a and study the following OP system 


I 


(TTfelTTfe/) = /“e “^A“7rfe(A)7rfe/(A)dA = hkdk,k> 


(3.1) 


7rfc(A) — A^ + ... , 

such that the norms hk’s are defined by 

hk = (TTfelTTfc). (3.2) 

As we show here, some useful relations can be obtained by varying a. Eventually, 
we will of course set a = 1 (2.2). The first polynomials can be computed from (3.1) 
to obtain 


7ro(A) = 1, 


7ri(A) = A - 


\ d 

a 


e 

r(l + a, to) 


(3.3) 

(3.4) 


where T{y^x) = y'^ ^dy is the incomplete gamma function. Obviously, the 

expression of the OPs becomes more and more complicated as k grows. The 
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polynomials tt^ being OPs, they satisfy a three-term recursion relation, which can 
be obtained as follows. As Att^ is a polynomial of degree fc -I-1, we can expand it on 
the basis of these OPs. Because {'Kk- 2 -i\^'^k) = 0 if i > 0, we can write the three 
term recurrence relation |44| : 


X'Kk '^k+1 “f 4“ ^k^k—1^ 


(3.5) 


where, by definition 


Sk^k — ('^/c I AtT/,;) , 

Rk^k—l — 11 Att^) . 


From (3.31 and (3.41, we can compute the first terms 

r(l -I- a, at) 


ho = /“ e-“^A“dA = 


-yfl + l 


1-1-a e 

So = -da log ho = -h 


a r(l -I- a, to) 


Ro — 0, 
Uo = 0. 


(3.6) 

(3.7) 


(3.8) 


3.1. Schlesinger equations 

In this section, we derive a couple of recursion relations called the Schlesinger 
equations that couple Rk and Sk- We first write 


Rkhk—l — 1 I AtT/c) — (AtT/,;—! [tt/,;) . 


(3.9) 


Therefore, using Eq. (3.5) with the substitution fc —>■ fc — 1, we have {XTTk-i\TTk) = 
i'^kl'^k) = hk and finally we obtain the standard relation 

hk 


Rk = 


hk -1 


(3.10) 


On the other hand, using the definition of the scalar product in (3.11, we have 

/ OO 

e““^A“ A7rfc(A)dA = -da{T^kWk) = -dahk, (3.11) 
from which we deduce the relation between Sk and hk 


Sk = -daloghk- 


(3.12) 


By combining Eq. (3.12) and Eq. (3.10), we obtain straightforwardly 
Sk+l - S'/C = -da log = - 


da log Rk+I- 


(3.13) 
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We now study the coefficient Cfc of the term of degree /c — 1 in the polynomial 

(3.14) 


[see Eq. (2.21]: 


'^/c(A) — + CfeA^ ^ + ... 

Taking the derivative of this equation with respect to (w.r.t.) a we obtain 

9Q7rfe(A) = 5aCfeA''“^ + ... . (3.15) 

Multiplying both sides of Eq. ( 3.15| ) by A and projecting on yields 

daQkhk = ■ (3.16) 


Besides, by looking at the term of degree k, i.e. oc A^, in Eq. (3.5) we find a 
recursion relation between Sk and 


Ck — Cfc+1 “t” ) 


which can be solved for C.k, using the initial condition (3.8 1 , to get 

k-l 

Ck = -Y.s^■ 


(3.17) 


(3.18) 


z=0 


Furthermore, by differentiating Eq. (3.111 w.r.t. to a we have 

/ OO nOO 

e““'^A“ A2 7r^(A) dA + 2 / e““'^A“ A i9ct7rfc(A)7rfc(A) dA 

= -(ATTfelATTfc) + 2(Aaa7rfc|7rfc) , (3.19) 

where, in the second line, we have simply used the definition of the scalar product 
in (2.2). Using the three-term recurrence relation (3.5) to rewrite the first term 

(3.20) 


and Eq. (3.16) to rewrite the second one, we have 

— hk-\-\ Sj^hk Rk^k — l “t” ‘^^aCkhk • 


We can also write 


— {daSk)kk Skidahk) • 


(3.21) 


By replacing the left hand side of Eq. (3.20) by Eq. dividing the resulting 

equation by and using Eqs. (3.10) and (3.12), we obtain a last recursion relation 
between Rk, Sk and 


Rk+l + Rk — ~9aSk + ^daCk ■ 


Combining Eq. (3.17) and (3.22), we have 


Rk-\-l Rk — ^aCk+l “t" ^aCk- 


(3.22) 

(3.23) 


And finally, with the initial condition i?o = Co = 0 in Eq. (3.8) we obtain 


^aCk — Rk • 


( 3 . 24 ) 
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Therefore substituting this relation in Eq. (3.221 we finally obtain a closed system 
of two coupled recursion relations between Sk and Rk 


^k+l ^a^ki 

Sk+1 lo§ Rk-\-l ; 


(3.25) 


where the second equation was previously obtained in Eq. (3.131. 

Our goal now is to relate a-derivatives to t-derivatives, i.e., find a relation 
between dahk and dthk- Differentiating Eq. (3.21 w.r.t. t, using that {dtTrk\T^k) = 0 
as dtTTk is a polynomial of degree fc — 1 (2.21, we have 

dthk = —e““*f“7r^(t). (3.26) 


We now start from the expression of d^hk given in Eq. (3.111 and use integration 
by parts to obtain 


/ OO 

clAir2(A)e-“*A“+‘ 

= [-ie-“^A“+‘4(A)]^; 


a 

1 

a Jt 

1 


dAe ((a + l)7r^(A) + 2A7rfc(A)9A7rfc(A)) 


= -{-tdthk + (a + l)hk + 2(7rfc|AaA7rfe)) , 
a 


(3.27) 


(3.28) 


where we have used Eq. (3.26). We can easily calculate the last scalar product 

(TTfelA^ATTfe) = (TTfcKfcA'' + ...)) = khk , (3.29) 

to obtain the desired relation between dahk and dthk 

adahk = tdthk — {2k + a + l)hk- 


(3.30) 


With this relation (3.30), it is then straightforward to relate daRk and daSk to 
dtRk and dtSk, using (3.10) and (3.12). Note that from now on, we set a = 1 

S'fe =-ti9t log/ife + (2A: + a + 1). (3.31) 


We finally obtain, using Eq. (3.25), the so-called Schlesinger equations 

Sk — Rk+i + Rk = tdtSk, 

2-Sk+i+Sk =tdt logRk+i- 


(3.32) 
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Note that, for a = 1, the initial condition (3.8) reads 
' ho = /“ e-^A“dA = r(l + a, t), 

So = —tdt log ho + a -\- 1 = Y{i+a,t) +0 + 1) 

Ro = 0, 

Uo = 0. 


(3.33) 


Using the Schlesinger equations (3.32) and this initial condition (3.33), we can 
compute step by step all the terms for arbitrary k. 


We end up this section by providing a useful relation between [see Eq. (3.141] 
d the cum 

t). One has 


and the cumulative distribution of the smallest eigenvalue Fiv(t) = Prob( min Ai > 

l<i<Af 


1^00 1^00 1^00 jyi 1 

'fiv(i) = / dAi / dA2... / dAjv Ujoint(Ai, A 2 , • • • = TT (3.34) 

Jt Jt Jt 

where the last equality is obtained by using the classical tricks of replacing the 
Vandermonde determinant by the determinant built from the OPs tt^’s (2.2) and 
then use the Cauchy-Binet formula mm- Using Eq. ( |3.18| ) with ( |3.31| ) as well as 

(3.35) 


(3.30) , we can write 


Cn = -N{N + a) + tdt log I /ifc j . 

\fc=o / 


And therefore, from the expression of E/v(t) given in Eq. (3.34), we have 

= -N{N + a) + tdt log (F^(t)) . (3.36) 

This expression thus provides a link between which is the first non trivial 
coefficient of the OP’s [see Eq. (3.14l[ and the cumulative distribution of the 
smallest eigenvalue. 


3.2. Laguerre-Freud equations 

In this section, we derive another set of recursion relations between the coefficients 
Rk’s and Sk’s, the so-called Laguerre-Freud equations, following the procedure used 
by Belmehdi and Ronveaux in |43j . To derive these equations, we start by searching 
two functions 'I' and <i>, which are polynomials in A, satisfying for any polynomial p 

(d/jp) = ($|p') , (3.37) 

where the polynomials $ and 4' may depend explicitly on the parameters t and 
a and p'{X) = d\p{X). Denoting by w{X) = e“'^A“ the weight associate to the 
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scalar product in Eq. (2.21 and using an integration by parts, this relation (3.371 
is satisfied provided that 

'^w + ($w)' = 0 and $(A = t)=0. (3.38) 

We find that the simplest non trivial solution to this equation is given by 
^'(A) = A^ — (2 + a + t)X + t{l + a), 


(3.39) 


$(A) =X‘^-tX. 


Then, using this property (3.371 for p = and for p = Tr^Tr^+i, we can write after 
expansion, using the three-term recursion relation (3.51 


^ 2 ,/c — (2 -|- a -|- t)Ii^k + t(l -l- a)/o,fe — 2 (J 2 ,fc — 

K2,k — (2 -|- a -|- t)Ki j^ + t(l -|- a)KQj. = L2,fe ~ 
where we have introduced the four family of integrals 

Im,k = (A-k^), 


Jni,k = (A'"|7rfe7r(,), 

^m,k ~ (A kfcll/c-t-l) 5 

^ Lm,k = (A'^KTTfc+iTTfc)'). 


(3.40) 


(3.41) 


These integrals are all calculated (for m = 0,1 and 2) in Appendix A Using these 
expressions together with Eq. (3.10|), we find the two relations 


k-l 


Rk+1 Rk T Sk(^Sk — ct — t — ‘2k — 2) -t- t(2k -t-1 -t- u) — 2 ^ 


i=0 

k 


(3.42) 


{Sk+i + Sk — ^ — a — t — 2k)Rk+i — 2 ^ ) Ri + ^ ) Sf — t ^ ) Si. 


i=0 


i=0 


We can rewrite these equations by subtracting rank k to the rank fc -|- 1 and find 
the two Laguerre-Freud recurrence equations (which are here of order 2): 

{ Rk +2 ~ Rk = Sk+i(2k -|- 4 -|- o -|- t — Sk+i) — Sk(2k + a + t — Sk) — 2t, 

5'fc+i(<S'fc+i — t) = Rk+i{2,k -|- 1 -f a -|-1 — Sfe+i — Sk) ^ ^ 

~ Rk+2{^k -|- 5 -f a -|-1 — 5^+2 — 5'fc+i). 

As we show below these two sets of equations ( 3.32[ ) and ( 3.43[ ) allow us to (i) derive 
the connection to the Painleve equation and (ii) perform the asymptotic analysis 
of these coefficients for large N. 


























13 


4. Painleve V equation for finite TV 

In this section, we proceed to the derivation of the Painleve V equation, following 
the method of Ref. [31]. First, it is useful to introduce the quantities |15] 

Ok — + 1 + a — Sk^ 

(4.1) 

— TTfc Cfc- 

Manipulating the Laguerre-Freud equations ( 3.43| ) we can prove (see Appendix B[ ) 
~ OkiOk — « — 2fc + t)ujk — Ok{kt(^k + a) + [Ok + t)Cfc) ~ 0 ; (4-2) 

which is a simple algebraic relation (and not a recursion relation) between the 
different variables ujk, Ok and (^k- Using the previous relations with the index k = N, 
we find from Eq. (3.361, 

Hn = tdt log(FAr(t)) = N{N + a)+ Cn- (4.3) 

Summing up the first Schlesinger equation (3.32[) from A: = 0tofc = TV — 1, we find 


N-l 


N-1 


N-l 




k=0 


k=0 


k=0 


which can be simplified using (3.181, (3.8l and (4.1l to yield 

— tdtCN = —Rn — Cn = ^N- 


On the other hand, using (4.3) and (4.5), we have 

= —wjv = Civ + Rn = H]S! — N{N + a) + Rn . 
Taking a derivative of this equation w.r.t. t, we find 

dtRN = t dfHN- 


From the second Schlesinger equation (3.32), we find 

tdtRN = Riv(2 — Sn — Sn-i) = Rn{0n — ^*iv-i) > 


where we have used (4.1). Using the relation derived in Appendix B in Eq. 
we can rewrite this recursion equation as a simple algebraic relation 

tdtRN = RnOn — TT^ • 

Un 


(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 


(4.9) 


Therefore, combining (4.7) and (4.9), we find 

t^dfllN = RnOn ~ 


On 


(4.10) 


On the other hand, from our Eq. (4.2) in which we inject the definition of ujn ( |4.1 
to eliminate Cat 


N(^N + a)t — (2TV + a)wjv ~ tRN = RnOn + 

On 


(4.11) 
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Summing and subtracting the last two equations (4.101 and (4.111, we obtain 


2^ = N{N + a)t - {2N + a)ujN - iRn - (4.12) 

On 

2Rn0n = N{N + a)t - {2N + a)wjv - tRn + (4.13) 


Multiplying these two equations (4.121 and (4.131 together, we find 


ARNtujf = {N{N + a)t - (27V + a)ujN - IRn^ - . (4.14) 


Finally, by using (4.6), we eliminate ujn and R^ (by expressing them in terms of 
Hn and dtHj^) and find the equation 


{tdfHN)^ = 4 {dtH^f (Hn - N{N + a)- tdtHN) + {{2N + a- t)dtHN + H^f , 

(4.15) 

which is a Painleve V equation in the Jimbo-Miwa-Okamoto a form |47I48] . Note 
that this equation coincides exactly with the equation first found by Tracy and 
Widom in [49]. 


5. Large TV asymptotic limit at the hard edge: Painleve III and 
first correction 

In this section, we study the behavior of the quantities Tijvj Rni Sn and in 


the hard edge limit (2.10), which, in the language of OPs, corresponds to a double 


scaling limit. Of course, as we are eventually interested in the study of the cumula¬ 
tive distribution of the smallest eigenvalue F)v(t), we could perform this asymptotic 


analysis directly on the Painleve V equation (4.15), as done by Tracy and Widom 


in their original study of GUE |T]. But it turns out to be much more convenient, 
especially to extract the 1/7V corrections, to perform this asymptotic analysis on 


the Schlesinger (3.32) and Laguerre-Freud (3.43) equations. 


To understand the structure of the coefficients Hn, Rn, Sn and Cat in this 


double scaling limit (2.10), it is useful to study their behavior for small t (keeping 


TV fixed). For t = 0, the OPs tt^’s in Eq. (2.2) can be expressed in terms of the 
generalized Laguerre polynomials 




r(a + A: + 1) 
k\ 


E 

2 = 0 


(-X)* 


i J r(a -I- T -|- 1) 


(5.1) 


Hence, thanks to standard properties of Laguerre polynomials, we easily obtain 
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these coefficients, for t = 0, using (|3.2[), (3.10) and (3.31) as: 


('^) lt=o 

= Li^\x)ki{-ir 

O 

II 

= T(k + a + l)fc!. 

o 

II 

= k{k + a), 

*^*4=0 

= ‘2k T n T 1, 

o 

II 

= —k{k + a). 


(5.2) 


We now take the index k = N. When t is closed to 0, we can use Eq. (3.26) and 
find 

dthN = -t“e-‘7rjv(t)2 = -rL^^(0)2iV!" + o(r) = + o(t“) . 

(5.3) 

Using (5.2), we can integrate and find the first correction 


Hn = r(7V + a + l)iV!- 


t“+i[r(fV + a + 1)]2 


+ o(t“+'). 


(5.4) 


r(a + 2)r(a+ 1) 

We are interested in the scaling behavior when N goes to infinity, t goes to 0 with 
X = Nt finite (2.10). In this double scaling regime, the expansion above reads 


1-1 


hN = r(iV + a + l)Af! 

This suggests the ansatz 

Hn = T{N + a+l)N\ 


N Vr(a + 2)r(a + l) 


+ o(x“+^)) + o( 1 




N 


where the function / has thus the small x expansion, read from Eq. (5.5) 

^^1 + 1 


/(a:) = - 


r(a + 2 )r(o + l) 


+ o(x“+i). 


(5.5) 


(5.6) 


(5.7) 


We can introduce the ansatz (5.6) in the relations (3.10) and (3.31) which give Rn 
and iStv in terms of the function /. We obtain 

Rj^ = N{N + a) + xf'{x) - f{x) + o (1), 


1 


Sn =2N + a+I - -^xf'iyx) + o ( — ) . 


1 


N 


(5.8) 


The first non trivial terms in Eq. (5.8) xf'{x) — f(x) for i?jv and —{1/N)xf'{x) 
for Sm are necessary to compute the leading order of the cumulative distribution, 
Fao{x). To compute the first 1/fV correction to the limiting distribution, we need 
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to expand i?jv and Sn in Eq. (5.8 1 to the next order in 1/A^. One actually expects 
the following expansion 


= N{N + a) + ^ nix)N-^ + o (N-^) , 


2 = 0 


(5.9) 


= 2Af + a + 1 + ^ Siix)N-'^ + o (TV^) 


i=l 


where the first term in this expansion are given in (5.8 1 , i.e., ro(x) = xf'{x) — f{x) 
and si(x) = —xf'{x). Besides, using the exact relation dthj^ = —t“e“‘7r7v(t)^ [see 
Eq. ^], one can show that ri{x) = o{x) as well as Si[x) = o{x) when a; —>■ 0. 
One can check, in principle, the validity of this expansion (5.111 order by order 
in powers of 1/fV by injecting it in the Schlesinger equations (3.32). Proving it 
rigorously to arbitrary order j is however a hard task. However, here, we only need 
this asymptotic expansion up to order 0{1/N'^), which can be obtained explicitly 
as follows. To compute the second correction, we truncate the expansion (5.111 up 
to the N~'^ order and write 


Rn = N{N + a) + xf{x) - f{x) + 


S — 2 -)- (j 1 — 


Xf'{x) , S2{x) , 


N 


iV2 


+ o(7V-2). 


(5.10) 


In the hard edge limit N —>■ oo, t —>■ 0, keeping x = Nt fixed, we can obtain the 
expansion of Rn+i and Sn+i at the same order by replacing by iV +1 and using 
{N + l)t = a:(l + 1/A^) in ( |5^ 

Rn +1 = (N + 1){N + 1 + a) + xf'{x) - f{x) + jf{ri{x) + x‘^f”{x)) 

+ ^{r 2 ix) + xr[{x) - ri{x) + \x'^f"{x) + \x^ f” {x)) + o{N-'^) 

Sn+i = 2Af + a + 3 - jjxf{x) + ■^{s 2 {x) - xf{x) - x'^f”{x)) + oiN-"^). 


By injecting these two expansions into the Schlesinger equation (3.32), we find 
(using tdt = xdx) at the first non trivial order (p{N~‘^) for the first Schlesinger 
equation and at 0{N~^) for the second) 


S2ix) - xs2ix) + ri{x) - xr[{x) = ^{x'^f"{x) + x^f"{x )), 


2s2(a;) — xS2{x) — xr[{x) = —ax^f"{x) + kx^f"'{x )), 


(5.12) 


which can be solved as follows. First, by subtracting the first equation of (5.12) to 
the second one, one obtains 


S2(x) = ri(x) - I a + - j x^f"(x) . 


( 5 . 13 ) 
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By injecting this relation (5.131 in the first equation of (5.12), one finds that ri 

(5.14) 


satisfies the following equation 


ri{x) - xr[{x) = -- (x'^f'ix) + x^f''{x)) , 
which can be solved, using that ri(x) = o{x) as x 0, yielding 

n(x) = ^x^r(x). 


Consequently, from Eq. (5.13) one has 


S2(a;) = ■ 


(5.15) 


(5.16) 


We finally find the two first terms of the expansion of Rn and Sn, for large N, as 


a 


= N{N + a) + xfix) - fix) + —x^f'ix) + o — , 


2N 

Sn =‘2.N + a + l- ^xf'ix) - °^^x^f''ix) + o 


1 


N 

1 


(5.17) 


From the expansion of Sn in Eq. (5.171, we compute C,n in the double scaling limit 
(2.10), using Eq. (3.17), up to the second non-trivial order for large N 


Cat = -{N + a)N + fix) + wr^xfix) + o ( — 


2N 


N 


(5.18) 


Finally, from the expression of Ctv we obtain the large N expansion of Elv(f) in the 
hard edge limit (2.10), using Eq. (3.36), which is given by 


xda: log (^Fn {^)) = fi^) + + ' 


(5.19) 


Using the initial condition JW(0) = 1 and (5.7) we can rewrite this equation as 




T7 = exp 


/(w) 


dw+^/(a:) + o( 4 


2N 


N 


Expanding Eq. (5.20) up to first order in 1/fV, one obtains finally 

where F^ix) is given by 

lim Fn (^) = Fooix ), Foo(a;) = exp ( [ ^^^du 
N —^oo \ jV / \t/0 


(5.20) 


(5.21) 


(5.22) 


Hence we easily obtain the functional form of the 1/iV correction as given in 
Eq. (5.21), fully consistent with the conjecture in (1.11) made in [3H]. However, 
at this stage, we still need to find the equation satisfied by the function /, which 
enters in the definition of F^oix) in Eq. (5.22). To obtain this equation, we expand 
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Cat in (5.8 1 and (5.18) to yield 


+0 


On and ujn in Eq. (4.1), which is easily done from the expansions of Rn, Sn and 

(5.23) 

UJN = -xfix) - ^{x^f'ix) + xf{x)) + O . 

Finally, by injecting these expansions ( 5.18[ ) and (5.23) in (4.2) with t = x/N, we 
obtain, by canceling the first term, of order 0{N~‘^), in Eq. (4.2) that / satisfies a 
Painleve III equation 

(.xf'f + 4/'(l + f){xf - /) = (a/)^ (5.24) 

with the small argument behavior in Eq. (5.7). This result coincides with the 
one obtained, by a quite different method, by Tracy and Widom |27| (note the 
correspondence (j(s) = —/(s/4), where a(s) is the notation used in |17]). Note 
that for integer values of a, f{x) can be written explicitly in terms Bessel functions 


P5] [see Eq. (C.l)]. These results in Eqs. ( 5.21|), (5.22) and (5.24) yield the results 
announced in the introduction in Eqs. ( |1.8 ), ( |1.9| ) and (1.11). Finally, in Appendix 
[C] we present a comparison between numerical simulations of Wishart matrices of 
size = 50 and the asymptotic formula in Eq. (5.21) describing the first 1/A^ 
correction. 


6. Large N asymptotic limit at the soft edge: Painleve II and first 
correction 

We now turn to the analysis of the PDF of the smallest eigenvalue Amin in the case 
where a ~ 0{N), and we set a = aN. In this case, the density of eigenvalues has 
a single support on [Nx-,Nx+] [see Eq. ( |l.4[ )[, with x± = {■\/l + a ± 1)^ with a 
soft edge at both extremities (see Fig. [^. Therefore, one expects that Amin will be 
close to Nx-, while its fluctuations, of order are governed by the Tracy- 

Widom distribution for (3 = 2. In the soft edge limit, the large N analysis of FN{t) 


in Eq. (2.1) is more conveniently done directly on the Painleve V equation (4.15) 
m- Following this route, one can indeed show [50] that, for large N 


7V1/3 

Amin = NX- -X ■ 


m 




where m is given by m 


(1 + a)i/® 


( 6 . 1 ) 


( 6 . 2 ) 


(^/^TA-1)4/3 ’ 

and where x is distributed according to the Tracy-Widom distribution F 2 , i-e- 
Pr[x < s] = ^ 2 ( 3 ) where J^ 2 (s) is given by [T| 


7^2(5) = exp ( — / (x — s)q^{x) dx 


( 6 . 3 ) 
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Here q{x) is the so-called Hatings-McLeod solution of the Painleve II equation 


q"{x) = xq{x) + 2q^{x) ,with q{x) ~ Ai(a;) , for a; —)■ oo 


(6.4) 


where Ai(x) is the Airy function. The result in Eq. (6.11 can be equivalently written 
as 


FN{t) = fo + o(l) 


(6.5) 


which implies 

HN(t) = tdtlogFNit) = -{mx-)hoix)N‘^^^ + o{N‘^^^) , x = m .(6.6) 

with hf){x) = fQ{x)/fo{x) and where fo = F 2 , and where the tilde refers to the 


soft edge scaling li mit. By injecting this form ( 6.61 into the Painleve V equation 

one fine 

Ah^ix) 


satisfied by Hjq{t) (4.15) one finds that ho satisfies the following equation 

2 _ _ “ 

— a;/ig(x) -I- ho{x) 


(^o(= 


(^o(= 


= 0 . 


(6.7) 


Using the Painleve II equation (6.4), one can indeed check that ho{x) = q'^{u)du 


is solution of this equation (6.7). Note that to check this, it is useful to use the 
identity 


ho{x)= / q'^{u)du = {q'{x) f - {q{x))‘^ - x{q{x)f 


( 6 . 8 ) 


What is the first correction to the limiting form in Eq. (6.5) ? Unfortunately, 


in the soft edge limit, it turns out the Schlesinger equations (3.32) do not allow 


to determine easily this first correction - while they were very helpful in the hard 
edge scaling limit. An alternative way to compute this first correction is to analyze 


directly the Painleve V equation in (4.15). By inspection of this equation (4.15), 


we conjecture that the first correction to Eq. ( |6.5[ ) is of the form 

Hnit) = —{mx-) (Jio{x)N'^^^ + hi{x)N^^^^ + . 


(6.9) 


By inserting this expansion (6.9) in Eq. (^4.15) we obtain that hi satisfies the 


following linear differential equation 


2hih'ii + 2{ho + /io(3^o - ‘2x))h'i + h'f^h'i = 0 , (6.10) 



amplitude A can not be determined from this analysis. One way to determine it 
would be to analyze the OP system ( |2.2[ ) in the limit when t is far from the left edge, 
i.e., for {Nx-—t) ^ (corresponding to the left large deviation tail of Amin El]), 
and then match this result with the typical regime, for \Nx- —1| ^ 0{N^/^). This 









































20 


program was carried out in detail in a similar albeit different context, involving 
discrete OPs in Ref. |3S] . Since we are interested in the first correction, this actually 
requires a very precise (and tedious) analysis of this regime {Nx-—t) ^ which 
goes beyond t he sco pe of the present paper. Hence our result for the first correction 
hi{x) in Eq. (6.101 does determine this function only up to a constant. We have 


not found any simple solution to this equation (6.101, which could indicate that the 


corrections to scaling in this case are actually more complicated, as found recently 
in the case of real Wishart matrices 1351. 


7. Conclusion 

To conclude, we have provided a direct computation of the cumulative distribution 
FN{t) of the smallest eigenvalue of complex Wishart random matrices dll] ), for 
arbitrary parameter a > 0. This was done by studying a set of semi-classical or¬ 


thogonal polynomials as defined in Eq. (2.21 for which we derived (i) the Schlesinger 


(2.8 1 and (ii) the Laguerre-Freud (2.9) equations. By combining these equations, we 
showed that Fj^ (t) can be expressed in terms of a solution of a Painleve V equation 
(4.15), thus recovering the result of Tracy and Widom |35] using a quite different 
method. In the large N limit, E)v(t), properly shifted and scaled, converges to 
a limiting distribution F^{x) which can be expressed in terms of a solution of a 


Painleve III equation (5.24) in the hard edge limit (corresponding to a = 0(1)) and 
of a Painleve II equation (6.7) in the soft edge limit (corresponding to a = 0{N)). 
Furthermore, we have computed explicitly the first correction to the limiting dis¬ 
tribution when N —>■ oo. In the hard edge case (1.11), we confirmed a conjecture 
by Edelman, Guionnet et Peche in Ref. [3S]. In this case, the first correction can 
be simply understood as a correction to the scale of the fluctuations of Aminj see 


Eq. (1.12). On the other hand, in the soft edge limit, we found that this correction 


is a solution of a second order linear differential equation with varying coefficients 


(6.101. Solving this equation remains a challenging open problem, which could sug¬ 


gest that the first correction does not correspond to a simple shift or rescaling of 


the scaling variable, as found for the hard edge (1.12). This could be reminiscent 


of the result found for real Wishart matrices m and certainly deserves further 
investigations. 
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Appendix A. Some useful integrals involving the orthogonal 
polynomials 

To compute the terms which enter the Laguerre-Freud equations, we need to com¬ 


pute several integrals containing the OPs 7rfe(x) (2.2). To perform these computa- 
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tions, we will follow the method developed by Belmehdi and Ronveaux |4H| . using 
Turan determinants. We introduce four types of integrals 

Im,k = (A™|7r2), 


Jm.k = (A™|7rfc7r^), 
Km.k = (A^lTTfcTTfe+i), 


(A.l) 


, Am.fc = (A™|(7rfe+i7rfe)') , 

for m = 0,1 and 2, while k is an arbitrary integer. Using the recurrence relation 
(3.5), we can calculate the integrals Im,kS and the Km,kS: 

' Io,k = (Ikfc) = hk, 

h,k = {\H) = Skhk, 
l2,k = (A^k^) = {Rk+1 + + Rk)hk, 

(A.2) 

(1-kfc^fc+i) ^5 

Rl.k (Akfc^fc + l) + 

. K2,n = (A^kfc^fc+l) = i^k+l + <5'fc)/l/c+l- 
The terms Jo,fe) Ji,k and Lp.fc can easily be calculated using that tt^’s are monic 
OPs: 


' Jo,k = (ikfeO = 0. 

Ji,k = (Akfcirk = 


(A.3) 


, Lo,/c = (l|(7rfc+i7rfe)') = {k + l)hk. 

The three last quantities, J 2 ,fe,Ti.fc and L 2 ,fc, require more work. For this purpose, 
we introduce a new object: the Turan determinant ^ defined as: 

'3^k+l — '^k+2'^k '^k+1' 


(A.4) 


Using the three term recurrence (3.5), we find a recurrence for the Turan determi¬ 
nant 

^k+i — (Att^^i Rk+i'^k)'^k Hfc-t-i (Att/c SkTTk Rfcirfc—i),(A.5) 

= RkT^k+lT^k-l + {Sk — Sk+l)TTk+l'^k — Rk+lT^k^ (A.6) 

= Rk.% + Fk, (A.7) 

where we have introduced the auxiliary quantity Fk = {Sk — Sk+i)T^k+iT^k + {Rk — 

Rk+iW- 
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Finally, we can write, using = hk/hk-i and the initial condition = Fq 

k 

■%+i = hk ^ (A.8) 

- /-» ^ i 

i—O 

This equation allows us to rewrite 

k 

(A™|^fe+i) =hkY. - S^+i)Lm,^ + 2 (i?, - R,+i)Jm,i ), (A. 9 ) 

i=0 

where the prime denotes a derivative w.r.t. x. On the other hand, by using the 
definition of the Turan determinant we find 


= (A'"|(7rfc7rfe+2)') - 2J„,fc+i. (A.IO) 


We obtain finally 

k 

(A™|(7rfe7rfe+2)^) = 2Jm,k+l + hf; ^ {{Si — Si+i)Lm,i + {Ri — ^i+l)2Jmy) . 

i=0 

(A.ll) 

We now have all the ingredients needed to calculate the last three integrals. 

Using the three term recurrence and the orthogonality condition, we can express 
Ll^k 3-8 


Li.k — (A|7rj,_|_]^7rfe) — Sk{TTk\T^k+i) F Rk{'n'k-i\'^k+i)i 
= Sk{k + l)hk + Rk{Trk-i\TTi+i). 


(A. 12) 
(A. 13) 


From Eq. (A.ll) for to = 0 , we can express the last scalar product in terms of 
integrals that we have computed before. Using Eq. (A.3l, we find 

k— 1 k 


Li^k — Sk{k + l)hk + Rkhk-i ^{Si — Si^i){i + l)hi — S'i.(A.14) 

i—0 i—0 

Using the three term recurrence and the orthogonality condition, we can express 
J 2 M as 


J2,k = (A^lTTfeTTfe) = Sk{\\nkn'^.) + i?fe(A|(7rfe_i7rfc)') 

= •S'fc (A.15) 


Finally, we have from Eqs. (3.10), (A.3) and (A.14), 


k-l 


J2,k = hk{kSk + ^ Si). (A.16) 

i^O 

Similarly, using again the three term recurrence and the orthogonality of the family, 
we can express L 2 ,fe as 


L2,k = (A^KTIfe + lTTfe)') 

= (A^lTTfe+lTT^) + (A|7rfc+1 T^k+i) + >5'fe(A|7rfe7r^_|_j) + RhiXlirk-Wk+i) 

= khk+i + Ji,fc+i + SkLi^k + Afc(A|(7rfe_i7r/j+i)'). (A.17) 
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Using Eq. (A.111 for m = 1, we can express the last scalar product in terms of 


integrals that we have computed before. 


L2,k — (2fc + l)/ife+i + Skh/, 'y ^ Si + Rk‘2'khj; 


i=0 


k-l 


+ hk 


i=0 


(5, - 5,+i) ^ Sj + 2iiR, - i?,+i) 
j=o 


We can simplify the telescoping sums and find 


L2,k — (2fc + l)/ife_|_i + /ife j 2 S': 


i=i 


1=0 


(A.18) 


(A.19) 


Equipped with the computations of these integrals, we can now derive the Laguerre- 
Freud equations. 


Appendix B. Details related to the Laguerre-Preud equations 

The derivation presented here follows the one of Ref. [TB] where slightly differ¬ 
ent OPs were considered. To derive the Laguerre-Freud equations, we start from 


equation (3.421 : 


k-l 

Rk+1 + Rk + Sk{Sk — o, — t — 2fc — 2) -|- t{2k -|- 1 -|- a) = 2 ^ Sj, 

1=1 

k k k 

{Sk+i +Sk — ^ — a — t — 2k)Rk+i = 2 Rj -|- — t Sj. 

1=1 1=0 1=0 

It is useful to introduce the following variables 

9k = 2 k + 1 + a — Sk, 

^k — Rk Ck- 


(B.l) 


(B.2) 

(B.3) 

(B.4) 


{X + t — Sk){X + uJk+1 — ^k) = ^fc+iRfc+i ~ 9k-iRk + A(A + t — Sk + ujk+i — ^k) 

= 9k+iRk+i — 9k-iRk + A(A -l-1 -|- Rk — Rk+i) 

= i9k+i — X)Rk+i — (9k-i — X)Rk + A(A -I- 1) . (B.5) 


In terms of these variables, we can rewrite the first equation of (B.l I as 

Wfe+i + ojk = (t — Sk)dk, 

and the second equation of ( |B.1[ ) as 

(t SkXjiy^k — 9k—iRk 9k+iRk+i- 

Similarly, Eq. (|B.4|) can be written as 
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We also rewrite Eq. (B.3l as 

(—+ X{2k + a — t) + 1 + ojk+i + Wfc) = (A +1 — Sk){dk — A). (B.6) 


Finally, by multiplying together Eqs. (B.5l and (B.6l, we have, with the notation 
^fc(A) = ~^ + f + a — t) + LOk, 

(f^fe+i(A) — ^2fc(A))(nfe+i(A) + r2fc(A)) = {9k — X){{9k+i — X)Rk+i 

- {0k-i-X)Rk + X{X + t)). (B.7) 


For A = 0, we find 

^k + l ~ Rk+l9k+l9k Rk9k9k — 1- (^■^) 

Using the initial condition Rq = loq = 0, we can solve this equation and find 

Wfe = Rk9k9k-i ■ (B.9) 


For A = —t, we find 

^k+l{—t)'^ — ^k{ — t)^ = {(^k + t)i9k+l + t)Rk+l ~ (fi'fe-l + i){dk + t)Rk- (B.io) 

Again, using the initial condition i?o = Wq = 0> we can solve this equation and find 
A^fc(—= {dk-i + t){9k + t)Rk + ^2 ' (B-11) 

Finally, using the explicit expression of VLk{—t), this equation can also be written 
as 


(wfc - kt){LOk - (fc + a)t) = Rk{9k + t){6k-i + t) . 


(B.12) 


Thanks to Eq. (|B.9|, we can substitute 0k-i in (|B.12|). Using the definition of uJk 

k 

(B.13) 


in (B.2|, Eq. (B.121 yields a quadratic equation for ojk in terms of 9k and (^k 

^k ~ (^k{9k — CL — 2k + t)ujk — 9k {kt{k + a) + {9k + t)Cfc) “ 0 ; 


which yields the equation given in the text in Eq. (|4.2|. 


Appendix C. Numerical simulations 


In this section, we present a numerical check of our formula for the first finite N 
corrections to the limiting density of Amin in the hard edge case, i.e., keeping a 
finite while A —)■ oo. Our simulations have been carried out for integer values 
of the parameter a for which the function f{x) in Eq. (5.221 can be computed 
explicitly in terms of Bessel I functions |2j^| : 


f{x) = -X 


det {J^j — k+ 2 {‘^“\/^')]l<j,k<a 
det {Ij — k{2^^/x'){'i-^j^k<a 


a€N* 


(C.l) 


and f{x) = —X for a = 0. In Fig. ([^ we show a plot of \ogFtq{x/N) obtained by 
sampling 3.10^ independent Wishart matrices of size A = 50 (dots) for different 
values of a = 0,1,2 and 3. The comparison with logFoo(x), evaluated exactly 
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Order 1 


Order 2 




0 



-5 


Figure 2. Comparison between our formula ( |5.20| l using (dashed curves) and data obtained 

by sampling 3.10^ independent Wishart matrices of size = 50 (dots) for different values of 
a = 0,1, 2 and 3. The left panels presents a comparison between \og{FN(x/N)) and \og{Foo{x)) 
given in (5.22 i. The right panel illustrates the finite N correction by comparing N log ^ 
to ^f{x), see Eq. (C.2l. 


x/N) '' 
Foo(^) , 


from Eqs. (5.22) and (C.l) shows a good agreement between theory and numerical 
simulations. The main objective of these simulations is to test the formula for the 
first order correction, given in Eq. (5.211 which can also be written as 

'F^ix/N)' 


log 


Fooix) 




oil/N) 


(C.2) 


In the left panel of Fig. we show a plot of Vlog [Fm{x/N)/F ao(,x)] as a function 
of X and compare it to |/(a:), where f{x) is given in Eq. (C.ll. This compari¬ 
son illustrates the accuracy of the asymptotic formula characterizing the finite N 
correction to Foo{x). 


Bibliography 

[1] C. A. Tracy, H. Widom, Level-spacing distributions and the Airy kernel, Commun. 
Math. Phys., 159, 151 (1994). 

[2] C. A. Tracy, H. Widom, On the orthogonal and symplectic ensembles, Commun. 
Math. Phys., 177, 727 (1996). 

[3] For a review see S. N. Majumdar, Les Houches lecture notes on Complex Sys¬ 
tems (2006), ed. by J.-P. Bouchaud, M. Mfeard and J. Dalibard [arXiv: cond- 




























26 

mat/0701193 . Random matrices, the ulam problem, directed polymers & growth mod¬ 
els, and sequence matching 

[4] J. Baik, P. Deift, K. Johansson, On the distribution of the length of the longest 
increasing subsequence of random permutations, J. Am. Math. Soc. 12, 1119 (1999). 

[5] M. Kardar, G. Parisi, Y.C. Zhang. Dynamic scaling of growing interfaces, Phys. Rev. 
Lett., 56, 889 (1986). 

[6] K. Johansson. Shape fluctuations and random matrices, Commun. Math. Phys., 
209 ( 2 ), 437 (2000). 

[7] M. Prahofer, H. Spohn, Universal distributions for growth processes in 1-h 1 dimen¬ 
sions and random matrices, Phys. Rev. Lett. 84, 4882 (2000). 

[8] T. Sasamoto, H. Spohn, One-dimensional Kardar-Parisi-Zhang equation: an exact 
solution and its universality, Phys. Rev. Lett. 104 , 230602 (2010). 

[9] P. Calabrese, P. Le Doussal, A. Rosso, Free-energy distribution of the directed polymer 
at high temperature, Europhys. Lett. 90 , 20002 (2010). 

[10] V. Dotsenko, Bethe ansatz derivation of the Tracy-Widom distribution for one¬ 
dimensional directed polymers, Europhys. Lett. 90 , 20003 (2010). 

[11] G. Amir, 1. Corwin, J. Quastel, Probability distribution of the free energy of the 
continuum directed random polymer in l-h 1 dimensions. Comm. Pure and Appl. 
Math. 64 , 466 (2011). 

[12] S. N. Majumdar, S. K. Nechaev, Exact Asymptotic Results for the Bernoulli Matching 
Model of Sequence Alignment , Phys. Rev. E 72 , 020901(R) (2005). 

[13] C. Nadai, S. N. Majumdar, Nonintersecting Brownian interfaces and Wishart random 
matrices, Phys. Rev. E, 79 , 061117 (2009). 

[14] P. J. Forrester, S. N. Majumdar, G. Schehr, Non-intersecting Brownian walkers and 
Yang-Mills theory on the sphere, Nucl. Phys. B 844, 500 (2011). 

[15] K. Liechty, Nonintersecting Brownian motions on the half-line and discrete Gaussian 
orthogonal polynomials, J. Stat. Phys. 147, 582 (2012). 

[16] G. Biroli, J.-P. Bouchaud, M. Potters, On the top eigenvalue of heavy-tailed random 
matrices, Eur. Phys. Lett. 78 , 10001 (2007). 

[17] Eor a review see S. N. Majumdar, G. Schehr, Top eigenvalue of a random matrix: 
large deviations and third order phase transition, J. Stat. Mech. P01012 (2014) 

[18] M. L. Mehta. Random Matrices, Academic Press, 2kd Edition (1991). 

[19] P. J. Forrester. Log-gases and random matrices, Princeton University, Press, Prince¬ 
ton, NJ (2010). 

[20] C. A. Tracy, H. Widom. Nonintersecting brownian excursions, The Annals of Applied 
Probab., 17(3), 953-979 (2007). 

[21] G. Schehr, S. N. Majumdar, A. Comtet, J. Randon-Furling. Exact distribution of the 
maximal height of p vicious walkers, Phys. Rev. Lett., 101(15), 150601 (2008). 

[22] 1. M. Johnstone, On the distribution of the largest eigenvalue in principal components 
analysis, Ann. Statist. 29 , 295 (2001). 

[23] C. Nadai, S. N. Majumdar. A simple derivation of the Tracy-Widom distribution of 
the maximal eigenvalue of a Gaussian unitary random matrix, J. Stat. Mech., P04001 
( 2011 ). 

[24] Y. Chen, G. Pruessner, Orthogonal polynomials with discontinuous weights, J. Phys. 
A: Math. Gen. 38(12), L191 (2005). 

[25] P. J. Forrester, The spectrum edge of random matrix ensembles, Nucl. Phys. B 402(3), 
709 (1993). 

[26] J. J. M. Verbaarschot and 1. Zahed, Spectral density of the QGD Dirac operator near 
zero virtuality, Phys. Rev. Lett. 70 , 3852 (1993). 

[27] C. A. Tracy, H. Widom. Level spacing distributions and the Bessel kernel, Commun. 









27 


Math. Phys., 161(2), 289-309 (1994). 

[28] P. J. Forrester, T. D. Hughes. Complex Wishart matrices and conductance in meso¬ 
scopic systems: exact results, J. Math. Phys., 35(12), 6736 (1994). 

[29] G. Akemann, T. Guhr, M. Kieburg, R. Wegner, T. Wirtz, Completing the picture 
for the smallest eigenvalue of real Wishart matrices, Phys. Rev. Lett. 113 , 250201 
(2014). 

[30] G. Gyorgyi, N. R. Moloney, K. Ozogany, Z. Racz, M. Droz, Renormalization group 
theory for finite-size scaling in extreme statistics, Phys. Rev. E 81 , 041135 (2010). 

[31] G. Schehr, S. N. Majumdar, Universal Asymptotic Statistics of Maximal Relative 
Height in One-dimensional Solid-on-solid Models, Phys. Rev. E 73 , 056103 (2006). 

[32] N. El Karoui, A rate of convergence result for the largest eigenvalue of complex white 
Wishart matrices, Ann. Probab. 34, 2077 (2006). 

[33] I. M. Johnstone, Multivariate analysis and Jacobi ensembles: Largest eigenvalue, 
Tracy-Widom limits and rates of convergence, Ann. Statist. 36 , 2638 (2008). 

[34] I. M. Johnstone and Z. Ma, Fast approach to the Tracy-Widom law at the edge of 
COE and CUE, Ann. Appl. Probab. 5, 1962 (2012). 

[35] Z. Ma, Accuracy of the Tracy-Widom limits for the extreme eigenvalues in white 
Wishart matrices, Bernoulli 18, 322 (2012). 

[36] G. Schehr, S. N. Majumdar, A. Gomtet, P. J. Forrester, Reunion probability of N 
vicious walkers: typical and large fluctuations for large N, J. Stat. Phys. 150 , 491 
(2013). 

[37] J. Baik, R. Jenkins, Limiting distribution of maximal crossing and nesting of Pois- 
sonized random matchings, Ann. Probab. 41 , 4359 (2013). 

[38] A. Edelman, A. Guionnet, S. Peche. Beyond Universality in Random Matrix Theory, 
preprint arXiv, 1405.7590 (2014). 

[39] E. Basor, Y. Ghen, Painleve V and the distribution function of a discontinuous linear 
statistic in the Laguerre unitary ensembles, J. Phys. A: Math. Theor. 42 , 035203 
(2009). 

[40] G. Schehr, Talk given at FOGMT4, December 16, 2014, http://lptms.u-psud.fr/ 
gregory-schehr/files/2014/12/Talk_F0CM.pdf, (2014). 

[41] F. Bornemann, A note on the expansion of the smallest eigenvalue distribution of the 
LUE at the hard edge, preprint arXiv: 1504.00235 

[42] W. Hachem, A. Hardy, J. Najim, A Survey on the Eigenvalues Local Behavior of 
Large Complex Correlated Wishart Matrices, to appear (2015). 

[43] S. Belmehdi, A. Ronveaux. Laguerre-Freud’s Equations for the Recurrence Coeffi¬ 
cients of Semi-classical Orthogonal Polynomials, J. Approx. Theory, 76(3), 351-368 
(1994). 

[44] G. Szego. Orthogonal polynomials, American Math. Society, Providence R.I., 4th 
edition (1975). 

[45] A. P. Magnus. Painleve-type differential equations for the recurrence coefficients of 
semi-classical orthogonal polynomials, J. Comput. Appl. Math., 57(1), 215 (1995). 

[46] P. J. Forrester, N. S. Witte. The distribution of the first eigenvalue spacing at the 
hard edge of the Laguerre unitary ensemble, Kyushu Journ. Math 61 , 457 (2007). 

[47] M. Jimbo, T. Miwa, Monodromy preserving deformation of linear ordinary differential 
equations with rational coefficients: 11, Physica D 2, 407 (1981). 

[48] K. Okamoto, On the r-function of the Painleve equations, Physica F 2 , 525 (1982). 

[49] G. A. Tracy, H. Widom, Fredholm determinants, differential equations and matrix 
models, Gommun. Math. Phys. 163 , 33 (1994). 

[50] T. H. Baker, P. J. Forrester, P. A. Pearce. Random matrix ensembles with an effective 
extensive external charge, J. Phys. A, 31, 6087 (1998). 








28 


[51] E. Katzav, I.P. Castillo. Large deviations of the smallest eigenvalue of the Wishart- 
Laguerre ensemble, Phys. Rev. E, 82(4), 040104 (2010). 



